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(57) Abstract 

The present invention pro- 
vides a wind turbine with improved 
efficiency, reduced static as well as 
dynamic structural loads, reduced 
noise, and facilitating grid connec- 
tion with improved power quality. 
According to the invention these 
and other features are fulfilled by 
a wind turbine comprising a wind 
velocity measurement system for 
determination of air velocities in 
front of the wind turbine and hav- 
ing an output for provision of de- 
termined air velocity values. 
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1 

A WIND TURBINE WITH A WIND VELOCITY MEASUREMENT SYSTEM 

FIELD OF THE INVENTION 

5 The present invention relates to a wind turbine comprising a wind velocity 
measurement system for determination of air velocities in front of the wind 
turbine. 

BACKGROUND OF THE INVENTION 

10 

Wind turbines are becoming an important source of electric energy. The 
majority of wind turbines are connected to the electric grid. As the number 
of larger turbines increases the need for reliable local wind measurements 
for power curve determination and for advanced control becomes more 
15 important. 



Having a history of millennia and a century of obscurity following the 
introduction of fossil fuel engines, the wind power technology had its 
revival in the mid-seventies as a consequence of the oil crises. Various 

20 traditional applications such as water pumping and grain-grinding and small 
wind turbines for charging of batteries in remote areas have. in modernised 
form had their renaissance. However, a large fraction of the global energy 
supply is now delivered in the form of electricity and thus the major 
effort in wind power technology was placed in the development of 

25 electricity generating machines. Therefore, nearly the entire manufactured 
wind turbine capacity is now so-called grid connected wind turbines. 

The modern "stock- average" wind turbine has an aerodynamic rotor with two 
or - most frequently - three blades. The rotational speed of the rotor is 

30 one-half to one revolution per second; through a gearbox the speed is 

increased to a level suitable to an electrical generator which is typically 
of the induction type: the rotor of the (electric) generator is magnetised 
by the grid and once the connection -to the grid has been established the 
rotational speed is kept constant within 1-2% over the range of operational 

35 wind speeds, 3-25 m/s . The nacelle (machine cabin) with the gearbox and 
electric generator is mounted on a tubular or lattice tower through a 
motor-driven yaw system which by means of a wind direction sensor keeps the 
wind turbine facing the wind. When the blades are fixed to the rotor hub 
the machine is said to be stall-regulated because the blade profile angle 

4 0 to the wind has been selected so that the blades stall and loose power when 
the rated power of the electric generator is reached. If the blades can be 
pitched the machine is said to be pitch-regulated: when rated power of the 
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electric generator is reached, the blades are pitched to limit the 
aerodynamic power. 

Wind turbines operating with variable speed are also known. 

5 

The grid connected wind turbine is equipped with a control system providing 
smooth cut-in and cut-out of the grid, correct yaw positioning of the 
rotor, and pitch setting of the blades, stopping of the machine if voltage 
is too high or too low or if loads on the phases are not balanced, stopping 

10 of the machine if unacceptable vibrations are experienced or the 

temperature of gearbox or generator is too high etc. A number of sensors 
are connected to the control system; sensors monitoring the air flow may be 
a wind vane and a cup anemometer mounted on top of the nacelle - just 
behind the rotor and the flow is therefore rather disturbed compared to 

15 free-flow conditions. 

The grid-connected wind turbines built and sold commercially in the 
seventies had installed capacities of 20 to 50 kW and aerodynamic rotor 
diameters of 5 to 15 m. In 1997, the typical commercially available 
20 machines have installed capacities of 500 to 800 kW and rotor diameters of 
40 to 50 m. The next generation of machines presently being tested are 1 to 
1.5 MW. This urge to increase size probably stems from belief in w economy 
of scale", i.e. the larger the cheaper, and from the fact that it is easier 
to fit into the landscape one 1 MW machine than ten 100 kW machines. 



Cost of energy from wind power plants depends foremost on the wind energy 
resources, i.e. the kinetic energy flowing through the area covered by the 
wind turbine's rotor (per sec): 



30 where A* is the area swept by the rotor, p is the air density, u is the wind 
speed and f is the frequency distribution of wind speed. The -consequence of 
the available energy being a function of the cubed wind speed, is that if 
the annual mean wind speed at one site is 20% higher than at another site 
the available energy is |(1.2) 3 -1) = 70% higher. Thus, moving from a site 

35 with an annual mean wind speed of 5 m/s to a site with 6 m/s the cost of 
energy (potentially) drops 1/1.7 = 58% thus, illustrating the vast 
importance of the wind climate to the economy of a wind power plant. 

The other factors are cost of the wind turbines, operation and maintenance 
40 and the efficiency of the wind turbines. Over 20 years, improved 
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aerodynamic and structural design tools and production methods have reduced 
cost of energy - disregarding the wind climate - by 70%. By 1997 the cost 
of energy at European sites with high wind speeds was 0.03-0.04 USD/kWh, 
and at "medium- good" wind sites 0.05-0.07 USD/kWh. Reduction in cost of 
energy becomes increasingly difficult. However, it is expected that the 
next 30 years will bring a further reduction of 20-30%. 

BRIEF DESCRIPTION OF THE DRAWING 

Figure 1 shows a conceptual lay-out for a wind turbine with control based 
on upstream wind measurements, 

Figure 2 shows a normalised output power, Cp, (lower traces) for a wind 
turbine operating at fixed and variable rotational speed, respectively 
(upper traces ) , 

Figure 3 shows the additional energy due to turbulence, 
Figure 4 shows a basic closed-loop control system, 

Figure 5 shows a contour plot of shaft torque as a function of wind speed 
and pitch angle (WD34 °35 rpm) , 

Figure 6 shows an overview of a wind turbine from a control point of view, 

Figure 7 shows a nomalized output curve, Cp, curve for WD34 for a pitch 
angle of 1 deg, 

Figure 8 shows an open loop estimation of wind speed, 

Figure 9 shows a closed loop estimation of wind speed, 

Figure 10 is a blocked schematic of a laser anemometer system, 

Figure 11 shows the number of photons/sec scattered by particles in a 
single spatial mode versus the optical wavelength assuming that 1 watt of 
light power is transmitted in a counter propagating mode otherwise 
identical to the light collection mode. The parameter v is a characteristic 
parameter of the Junge particle size distribution. A large v implies a 
relatively large amount of small particles. Note that the number of photons 
per second peaks at around 10 um almost independently of v. The wavelength 
dependence of the refractive index has not been incorporated, 
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Figure 12 shows a layout of the optics of a back-scattering reference beam 
configuration. The polarising beam-splitter in conjunction with the wave 
plates ensures efficient utilisation of the available laser power, 

5 

Figure 13 shows a layout of an autodyne or self -mixing configuration. The 
scattered light propagating in the same mode as the transmitted light - 
except for the direction of propagation - enters the laser and perturbs the 
internal power level. This is detected by the photodetector at the rear end 
10 of the laser. It is noted that a it is very important to avoid reflections 
that would enter the laser in a mode that could be amplified in the laser, 

Figure 14 shows an autodyne system with a tapered waveguide laser, 

15 Figure 15 shows a beam expander based on reflective optics. The location of 
the measuring volume in a plane perpendicular to the initial optical axis 
can be changed by tilting the largest of the two mirrors. The distance to 
the measuring volume can be changed by a displacement of the large mirror. 
It is important that it is the large mirror that is moved. 

20 

Figure 16 shows a signal processing layout. The frequency locked loop 
(inside the dashed line) may be implemented in the form of an application 
specific integrated circuit. The microprocessor sets the range as well as 
the loop speed and will also control the initiation of the lock. The input 
25 filter is envisaged to be based on traditional analogue circuitry. It 
serves to eliminate the dc signal from the detector and to confine the 
noise bandwidth, 

Figure 17 shows the beam parameters: d t is the diameter of the 
30 transmitter/ receiver aperture, d c is the focal beam diameter, L is the 

distance to from the anemometer to the centre of the measuring volume and 
2Az is the length of the measuring volume, 

Figure 18 shows a simulated detector signal assuming 10 % turbulence, 
35 The third line concerns the development of advanced control strategies and 
systems that will improve efficiency, reduce static as well as dynamic 
structural loads, reduce noise, and facilitate the grid connection in terms 
of power quality. Also the application of power electronics, allowing for 
variable-speed operation, will probably be common in the future, 

40 

Figure 19 shows characteristic parameters in conjunction with velocity 
measurement from a wind turbine, 



SUBSTITUTE SHEET (RULE 26) 



WO 98/42980 



5 



PCT/DK98/00125 



Figure 20 shows plots of H R (k) for various weighting functions across the 
rotor, 

5 Figure 21 shows a plot of the correlation function a(Xii), 

Figure 22 shows a plot of the correlation function a°(Xoti 2 ), and 

Figure 23 shows a plot of the error function erf (x) . 

10 

SUMMARY OF THE INVENTION 



It is an object of the invention to provide a wind turbine with improved 
efficiency, reduced static as well as dynamic structural loads, reduced 
15 noise, and facilitating grid conn ection with improved power quality. 

According to the present invention the above and other objects are 
fulfilled by a wind turbine comprising a wind velocity measurement system 
for determination of air velocities in front of the wind turbine and having 
20 an output for provision of determined air velocity values. 

The wind velocity measurement system may comprise a laser anemometry system 
having a laser for emission of a laser light beam for illumination of 
particles in the air, a detector for detection of light emitted from 
25 particles that have interacted with the illuminating laser light, and a 

signal processor connected to the detector for generating a velocity signal 
corresponding to the velocity of the particles based on the detected light. 

According to a preferred embodiment of the invention, a wind turbine is 
30 provided comprising adjustable pitch blades, a pitch controller for setting 
blade pitch and having an input, the blade pitch being set according to a 
signal provided at the input, and a control system for controlling the 
output power generated by the wind turbine by control of blade pitch. The 
control system may have a control output that is connected to the input of 
35 the pitch controller, a velocity input that is connected to the output of 
the velocity measurement system, and a pitch processor for calculation of 
the desired blade pitch based on determined air velocity values and 
generation of a corresponding pitch control signal that is provided at the 
control output. 

40 

It is an important advantage of the invention that the determination of air 
velocities in front of the wind turbine compensating for delay in blade 
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pitch setting so that the actual blade pitch setting corresponds to the 
actual wind velocity at the blades during operation of the wind turbine. 

According to another embodiment of the invention, a method of wind turbine 
5 power curve determination is provided comprising the steps of mounting on 
the wind turbine a wind velocity measurement system for determination of 
air velocities in front of the wind turbine, connecting a power measurement 
system for determination of the output power generated by the wind turbine 
between a power generator of the wind turbine and a load receiving the 
10 power generated by the wind turbine, and recording corresponding values of 
wind velocity from the wind velocity measurement system and output power 
from the power measurement system. 

According to still another embodiment of the invention, a method of 
15 controlling output power from a wind turbine is provided comprising the^ 
steps of mounting on the wind turbine a wind velocity measurement system 
for determination of air velocities in front of the wind turbine and 
providing determined velocity values to a power control system of the wind 
turbine. 

20 

Optimisation of wind turbine efficiency - and in turn measurement of the 
efficiency - includes accurate measurement of wind speed in the immediate 
vicinity of the wind turbine. A laser anemometer mounted on the wind 
turbine could serve both purposes possibly with better accuracy than the 
25 presently used cup anemometer, see figure 1. 

Wind turbines may have pitchable blades and operating with variable 
rotational speed of the rotor providing optimal energy extraction and 
improvement of power quality. Under given flow conditions - mean wind 

30 speed, turbulence shear etc. - there is an optimal setting of rotor speed 
and blade pitch. However, the response time of the wind turbine is several 
seconds and thus, a short-term prognoses of the flow characteristics is 
needed. The laser anemometer may potentially deliver such forecast. The 
characteristics of a wind turbine operating at a constant and variable 

35 rotational speed, respectively, is illustrated in figure 2. The term C p is 
a normalised power output of the wind turbine, defined by the following 
expression: 

where p(u) is the so-called power curve of the wind turbine: the power out 
40 as a function of wind speed. The efficiency of fixed-speed wind turbine is 
a function of wind speed and is maximum only in a narrow wind speed band. 
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Variable-speed operation will expand the range of high (maximum) 
efficiency, and the laser anemometer facilitates optimal control of speed 
with regard to efficiency - and potentially also for structural load and 
noise reduction. 

5 

Further, the laser anemometer facilitates (retrospective) assessment of 
production capability of wind power plants, being important to the 
credibility of the wind power technology in the financial world. The 
measurement of wind turbine efficiency involves measurement of wind speed 

10 at approximately 50 m above the ground and 100-200 m upstream, where the 

flow is undisturbed of the machine itself - figure 1. The laser anemometer 
may measure the wind speed without having to erect a separate tower; in 
addition, the laser anemometer may provide a spatial average of wind speed 
and turbulence accurately representing forces on the wind turbine. This 

15 would contribute to an improvement of the accuracy of performance 
assessment of wind power plants . 

Below, characteristics are outlined of the boundary layer flow field. Such 
characteristics are assumed important in determining the performance of a 
20 wind turbine and in selecting a strategy (for maximum energy extraction) 
for the rotational speed of rotor and/or blade pitch angel. 

Power output from a wind turbine may be predicted based on knowledge of the 
wind speed frequency distribution at the wind turbine site and of the so- 

25 called power curve of the wind turbine. The frequency distribution shall be 
given for the wind turbine's hub height, and the power curve is measured 
with an anemometer elevated to hub height upstream of the wind turbine. In 
both cases 10 min averages of wind speed (and power) are applied. As the 
same averaging time is used for both power curve and frequency 

30 distribution, the end-result is independent of averaging time. However, in 
the real world turbulence, flow shear, air density, etc, may be different 
for the site where the power curve was measured and the site on which the 
wind turbine is to be positioned. Experience indicates that in fact these 
parameters are important and uncertainties are introduced when disregarding 

35 them.. 

As stated, when measuring power curves turbulence is disregarded as being 
of lesser importance. However, in stationary flow with mean U and (for 
illustration) isotropic turbulence (the flow has identical characteristics 
4 0 in all directions, thus also the same fluctuations in the three directions: 
au = av = ow) the kinetic energy in 1 m 3 of air is 
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E = -p(U+uf+ v 2 +w 2 =-pU 2 + -pat (3) 

where overbar denotes average (note that it is assumed that the air density 
is constant) ; the energy flux (power) through an area of 1 m 2 perpendicular 
5 to the mean wind direction is 



P=(U + u)E = Vi p [(U + u f +v 2 + \v 2 ](U + u)*Vi pU 3 0 + 5 ilk <*> 

where I u = a u /U is the turbulence intensity. Thus, for turbulence 
10 intensities 10-20% the contribution to "per-m 2 -energy-f lux" from turbulence 
is 5-20%. More realistically, in the atmosphere, the coefficient to I u is 
of the order 3-4 rather than 5 since au is larger than the other two 
components.- Thus, for turbulence intensities 10-20% turbulence .increases 
the energy flux per m 2 with 5-15% additional to what is calculated from the 
15 mean wind speed alone. Disregarding lateral and vertical components of 
turbulence the energy flux becomes 



P = Vi p(U + u f (U + u)*V2 pU* 0 + 3 il), (5) 

20 Here, the contribution to energy flux from turbulence of 10-20% is 3-12%. 
Performance data indicate that power output may vary up to 8 % for 
variation in turbulence intensity from under 5% to over 15%. In performance 
standards, it is presently not taken into account that the energy flux 
depends on turbulence intensity. 

25 

Thus, improvement of the accuracy of performance testing may amongst other 
things be improved by also measuring the turbulence that the machine is 
exposed to. To do that a suitably spatial/time average of turbulence must 
be defined and measured. The laser anemometer facilitates contribute 
30 measurement of turbulence. 

Evidence points to higher power output when turbulence increases even for 
fixed-speed machines. If the rotor speed can be varied a larger fraction of 
the turbulent energy may be extracted, and it may be assumed that a 
35 continuous prediction of the instant wind speed - such as supplied by a 
laser anemometer - may help select the optimal rotor speed at any given 
time . 

Considering the potential improvement of energy extraction by variable- 
40 speed operation of the wind turbine, it is useful to investigate the 
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potential gain. The gain may be divided into what is gained by a relative 
slow regulation (time constant of the order of 10 min) and a faster 
regulation (time constant of a few seconds). 

Usually, climatology wind statistics is based on 10 minute average values 
of the wind speed. It is interesting how the averaging time influence the 
estimation of energy flux. Considering only the U-component of the wind 
field, the annual average per-m2 energy flux is 

P = F*(k)A 3 = V2pr(l + j)8.76A 3 



where A and k are the Weibull parameters and F E a function tabulated in 

Riso Report 1981 by Petersen et al. If the .turbulence intensity (I u = cr u /U) 

is zero, the energy flux is given by eq. 6 applying values of A and k 
determined from measured 10 minute averages of wind speed. If the 
15 turbulence intensity is not zero, the annual energy flux will be under- 
estimated. The difference in Weibull parameters between 10 minute averages 
and instantaneous values of wind speed has been evaluated. The relative 
differences, 5A/A and 5k/k, respectively, are for small values 

5k kil SA 1 n l ,3k 

20 — = — , a — ~-rV(l + -r)— (7) 

k 2(y/(l+2/k)-y/(l + l/k)) A k k'k 

where *F is the digamma function: 

, % d In Y(x) 
dx 



It is seen that 8A and 5k are functions of k and Iu. The energy flux 
neglected when using 10 minute average values is 

SP _ Fs(k + 8k) ■ A 3 (k + 8C) - F E (k) ; A? to 
P ~ F E (k)A 3 (k) 

For values k = 1.5 and 3, the eq. 9 is plotted in figure 3 together with 
the expression of eq. 5. 
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The points are the increase calculated by means of eq. 9 and the solid line 
is the simple eq. 5. It is seen that the energy flux increase is only 
slightly dependent on k and that the eq. 5 gives a good estimate not only 
in stationary flow but also in the atmospheric boundary layer, of the 
5 fraction of energy flux left out when using 10 minute average as a 
climatological base number. 

At realistic turbulence levels of 13 to 18% the energy flux ignored by 
using 10 minute statistics to calculate frequency distribution parameters 
10 is 5 to 10%. 

Since the turbulent kinetic energy usually is ignored, the non-constant 
turbulent kinetic energy appears as an intrinsic uncertainty, e.g. an 
uncertainty which is not eliminated by improvement of test equipment or 
15 duration of test. 



Thus, in highly turbulent air flows, the energy flux can not be accurately 
determined from the mean wind speed only. However, there are other 
intrinsic contributions to be considered: the wind turbine extract energy 

20 from an extended area (the swept area of the rotor), the single-point 

measurement at hub height may not sufficiently well represent the flux of 
energy experienced by the rotor. Thus, also horizontal and vertical wind 
shear and misalignment (yaw error) may cause errors in the power curve 
measurement. It has been found that shear and yaw error may cause 

25 variations in predicted annual energy output of more than 5%. 

The largest, present-day wind turbines have two or three blades of lengths 
up to 40-45m; the rotational speed of the blades is high, with blade tip 
speeds of the order 60 m/s. The blade profiles are similar to airplane wing 
30 profiles, the purpose of which is to reduce drag and in turn energy 
conversion losses. 

For the ideal rotor - e.i. a rotor with no drag losses and with infinite 
high rotational speed - operating in a non-turbulent flow, it can be shown 
35 that the maximum of the rotor efficiency defined as 

P 

C P = ^7— dO) 

' Kpu'Ar 

(where P is extracted power, A r the swept rotor area, p air density and u 
40 wind speed) is 
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(11) 



Further, it can be shown (still employing simplest possible theory) that 
the maximum efficiency is obtained when the flow speed behind the rotor is 
reduced to 1/3 of the upstream flow speed and to 2/3 at the rotor disc. 
Betz theory is inadequate when dealing with real rotors; however, it does 
point to the important fact that not only is the flow as measured by the 
presently used nacelle cup anemometer disturbed by the nacelle itself but 
the flow speed is inherently linked to the performance of the rotor: e.g 
the higher performance the lesser wind speed is measured on the nacelle 
top, implying that the value of this instrument for performance 
measurements is limited. 

Presently, the so-called blade element theory is the method applied for 
aerodynamic load -and power calculations, though with various auxiliary 
methods to take into account e.g. turbulence, stall behaviour of the blades 
and the induction or wake delay. The blade element theory quite accurately 
models power output and loads under stationary and quasi-stationary 
conditions when the flow around the airfoils are far from stall. The blade 
element theory does not account for dynamic flow behaviour. Therefore, to 
model the time-wise development in aerodynamic power, modifications have to 
be applied. 

Below, the conventional method for power performance measurements is 
reviewed and intrinsic uncertainties which may be reduced or eliminated by 
employment of nacelle mounted laser anemometry are pointed out. 

Since the IEA recommendation for power performance testing was issued in 
the late 1970' ies it has been found that actually unacceptable large 
uncertainties were attached to the tests : even under fairly ideal test 
conditions is was estimated that the uncertainty could reach a level of 10 
to 20%. 

Due to lack of a genuine commercial market this fact was not seriously 
acknowledged until a rapidly increasing commercialisation has hit the wind 
power industry in the resent few years. In reality, firm contracts between 
manufacturers and developers now implicitly call for uncertainty levels for 
performance testing as low as 1-2% in terms of standard uncertainty, which 
is unrealistic - even though the understanding of the problems has 



improved. 
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As of now a "Committee Draft" for power performance testing, IEC 96, has 
been prepared. The document has not yet obtained the status as an 
international standard but will probably do so in the near future. The 
document offers a procedure to be followed when measuring power curves, as 
5 well as a method for assessing uncertainties. The method and the proposed 
uncertainty analysis are outlined below. 

While the measured power curve may be used for other purposes, the IEC 
document is directed toward application to calculation of annual energy 
10 production (AEP) . The power curve is basically established by collecting 
pairs of wind speed and power measurements, though also additional 
measurements are made. 

Test site: shall be assessed for sources of flow distortion, mean terrain 
15 or other features causing differences in wind speed at the position of the 
wind turbine and the met towerT~At~ least topographical variations, other 
wind turbines, and obstacles shall be taken into account in the analysis. 

Meteorological tower: shall be placed between 2 and 4 rotor diameters (D) 
20 from the wind turbine; more than 2D to avoid stemming disturbances, and 
less than 4D to maintain good correlation. 2 . 5D is recommended. 
Measurements are acceptable for approx. 130° around the case where the wind 
is from met tower to wind turbine. 

25 Correction factor due to flow distortion: if the test terrain adheres to 
specified requirements no corrections need to be made and the uncertainty 
related to flow distortion shall be set to not less than 2-3% depending on 
separation between met tower and wind turbine. If the test terrain does not 
adhere to the specified requirements of an ideal test terrain site 

30 calibration shall be performed. If a computer model is applied the 

uncertainty shall be assumed to be 1/2 of the correction made to wind 
speed; if measurements used for site calibration the uncertainty is assumed 
to be 1/3 of the correction. 

35 Electric power: The sensor (s) to measure net electric power shall be of 
class 0.5 or better. 

Wind speed: Shall be measured with a cup anemometer mounted at the wind 
turbine's hub height H h "2.5% above the ground in the met tower. It is 
40 specified how to mount the instrument so as to avoid flow distortion from 

the met tower itself. The distance constant of the anemometer shall be less 
than 5 m/s, the instrument shall be calibrated and the uncertainty stated. 
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Wind direction: Shall be mounted at H h "10% in the met tower and the 
accuracy shall be better than 5%. The purpose of the instrument is to 
single out data, where the met tower is in the wake of the wind turbine. 

5 Air density: Should be (indirectly) measured at H h at the met tower, by 
measurement of pressure and temperature in accordance with ISO 2533. 

Precipitation: Shall be monitored to distinguish between wet and dry 
periods . 

10 

Wind turbine status: A parameter indicating whether the wind turbine is 
ready for operation or not. 

Data acquisition system: Sample frequency shall be at least 0.5 Hz. It 
15 shall be ensured that uncertainties related to the data acquisition system 
are negligibfe compared to other uncertainties. 

Measurement procedure: During the measurement the wind turbine shall be in 
operation and the above parameters shall be sampled at 0.5 Hz or faster. 
20 From the data mean values of 10 min periods of consecutive data sets, 

standard deviation, maximum and minimum shall be computed and stored. At 
least 180 hours of data shall be collected, and there must be at least 30 
min of data in each 4 m/s wind speed bin. 

25 Data correction: The data of the data base must be corrected (calibration 
errors, site calibration correction etc.) 

Data normalisation: Data shall be normalised to two reference air 
densities: one rounded off to the nearest 0.05 kg/m 3 and one to the ISO 
30 standard atmosphere 1.225 kg/m 3 ; the precise method is given in the IEC 
document. 

Power curve: The corrected and normalised data are binned in 2 m/s wind 
speed bins and averages of wind speed and power are formed in each bin. 
35 These sets of bin-averaged values of wind speed and power constitute the 
power curve. 

Annual energy production (AEP) : The AEP is estimated by applying the 
measured power curve to annual mean wind speed, assuming the lOmin average 
40 wind speeds to be Rayleigh distributed. The AEP is what in metrology is 
termed the measurand. 
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Knowledge - or at least awareness - of the component uncertainties and 
their mutual relations in an experiment is the key to improvement of the 
test. 

A guide to uncertainty analysis, ISO 1993, was issued recently, having all 
5 major standard organisations behind it. One of the main points of the guide 
is to disregard the concept of bias uncertainty and replace it with 
"correlated uncertainty". Doing so, all uncertainties may be treated 
equally, e.g. added in the same way, which is a major improvement. However, 
it somewhat complicates the presentation of the issue. 

10 

The true value (the concept of "true value 7 ' is also estranged in ISO 1993; 
it is used in the present disclosure anyway) of the measurand is never 
exactly known, and therefore determination of the uncertainty is a 
probabilistic problem, i.e. what is the probability that the true value is 

15 found inside the interval [u-a;u+a] , where p. is the measured value and a is 
the uncertainty measure? In some cases objective measures of the sought 
probability exist, i.e. when a stochastic variable is measured (with no 
instrumental uncertainty) and the measured values are independent and form 
a Normal Distribution with a mean u. and a standard deviation a, the 

20 probability of each measured value being in the interval [ji-a;u-kjj is 

approximately 68%, the mean and standard deviation being estimated from the 
measurements. In other cases the shape of the distribution and/or the 
standard deviation may not be explicitly known. 

25 If the distribution is not Normal, the probability attached to the standard 
deviation may be different. However, in this annex, it is assumed that all 
uncertainty components are sufficiently accurately characterised by a 
standard deviation. Further, it is trivial - though important - to note 
that the standard uncertainty for e.g. a random variable is not identical 

30 to standard deviation of the process, a. The uncertainty standard deviation 
is the estimated standard deviation of the final average {binwise or else- 
wise) of a component in contrast to the standard deviation, a, of the 
parameter. 

35 There are two types of uncertainties. It is assumed that they are only 
distinguished by the way their magnitudes are determined. 

The first type of uncertainty, category A, is characterised by firm 
knowledge of the uncertainty distribution, e.g. through the data base of 
40 the power performance test or tests made prior to the performance test. 
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The distribution of the measured quantity is assumed Normal and the 
standard uncertainty is estimated by the expression 



where a is the standard deviation of the measured process and N is the 
number of observations. Uncertainties of category A are in principle 
independent of the person who makes the uncertainty estimate, and the 
uncertainty of the estimate of the mean value can in principle be reduced 
to any small quantity by repeated measurements. 

The second type of uncertainty, Category B, is characterized by no or 
imprecise knowledge of the uncertainty distribution. Category B 
uncertainties is estimated by other means than measurements. The 
uncertainty type would typically be associated with "bias" uncertainty; 
however,, if a bias/one-sided error is identified the- -data- should be 
corrected. What is left after correction may or may not be biased; in any 
case the sign is not known. 

The (important part of the) uncertainties of instruments is often of 
Category B: little or no specific information may be available - or a 
statement of uncertainty limits, for instance f, U, is given in data sheets 
for the instrument. For instruments that are only specified by uncertainty 
limits, an uncertainty distribution as well as uncertainty standard 
deviation shall be estimated. Often, a rectangular distribution is chosen. 
Uncertainties of category B are often based on personal judgment, and 
consequently they are often dependent on the person who makes the 
uncertainty analysis. Uncertainties of category B can not be minimized by 
repeated measurements . 

The square of the standard uncertainty of annual energy production (AEP) 
can in its general form be expressed as 

N N M M 



a 



u = 



(12) 




(13) 



i=l ;'=/ k=l ;=/ 



where 



fi 



= relative occurence of wind speed in wind speed bin i, 
= standard uncertainty of component k in bin i, 
= sensitivity factor of component k in bin i, 



u k ,i 



Clc f i 



M 



= number of uncertainty components, 
= number of bins, 

= 8760, number of hours per year, 



N 
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Pk,ia,j = correlation coefficient between uncertainty component k in bin i 
and uncertainty component 1 in bin j . 

For all uncertainty components it is assumed in the IEC document that 
5 • different uncertainty components are uncorrelated, pk,i,i,j = 0 for k * 1. 

In addition, for type A uncertainties it is assumed that 
• Uncertainty components in different bins are uncorrelated, Pk,i,i,j - 0 
for i * j . 

10 

This implies that the square of the integrated type A uncertainty can be 
written as 

^=^SZ fV*A»l < 14 > 

1=1 k=l 



15 where M A is the number of uncertainty components of type A and A u 2 i,j< is the 
uncertainty of type A f component k in bin no. i. 

For category B uncertainties, the additional - possibly conservative 
assumption is that bin-wise uncertainties of the same component is fully 
20 correlated throughout all bins, pk r k,i,j = 1- 

This implies that the square of the integrated type B uncertainty can be 
written as 



25 «l=^Z(t F < 15 > 

*=i i=i 

where M B is the number of uncertainty components of type B and B u 2 i,k is the 
uncertainty of type A, component k in bin no. i. The motivation for the 
assumption of the type B uncertainties each being fully correlated over all 
30 bins is that a typical type B uncertainty is the uncertainty of an 

instrument's calibration constant, which presumably remains constant during 
the test. The total standard uncertainty on annual energy production is 



35 



U AEP = V^ +W * (16) 

In performance testing as in many other tests the type A uncertainties are 
less important since these are reduced any small number if the duration of 
the test is long enough, Eq. 12. However, the type B uncertainties remain 
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constant once test method, equipment, site etc. have been chosen. 
Therefore, improvement of the accuracy of performance testing primarily 
implies reduction of the type B uncertainties . 

Below, the uncertainties constituting important error sources are listed 
and improvements of the performance test provided by the laser anemometer 
are assessed. What presently appear to be the contributing uncertainties in 
a performance test is outlined in the table below. The list of uncertainty 
components may well be incomplete, however it reflects the present 
knowledge. The characteristics of main groups of uncertainties are outlined 
in the notes to the table. Typical ranges of uncertainties for each 
uncertainty component are listed in the table. 

The laser anemometer appears to be capable of providing improvement options 
to all uncertainty components, except "obstacles" : 



Instrumental (anemometer) : There are indications that calibration accuracy 
of (cup) anemometers in wind tunnels at present cannot be made better than 
1-2% due to blockage effect and lack of natural turbulence. The cup 



turbulence both U and V are non-zero, in the wind tunnel V=0. Therefore, 
the conditions of the wind tunnel are intrinsically different from the real 
world, which introduces an error. The laser anemometer is a "first- 
principle" instrument, which in principle does not need calibration, and it 
will basically measure only the U-component of wind speed, presumably the 
component experienced by the cosine-responding wind turbine propeller. 

A wind turbine mounted laser anemometer may facilitate experimental site 
calibration, both in terms of accuracy and in amount of data gathered 
within a certain period of time. To apply the laser technology for site 
calibration the instrument must be capable of measuring (at least) the 
axial wind speed at selected distances from the wind turbine. Also, with 
the availability of an operational laser anemometer, site calibration will 
not be a separate test but an integrated part of the performance test. 

Assessing the production capability of a wind farm the power curves are 
measured on a sample of wind turbines, typically one or two. With the 
achievable accuracy of the individual power curves, larger samples are 
needed. However, the per-unit present costs of power curves are high (due 
to met towers and instruments) . A low-cost laser anemometer would make it 
possible to performance test all units, and thus removing sample 
uncertainty completely. 



anemometer measures the horizontal speed, 




While in the natural 
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Amount of data: The laser anemometer may collect data from a larger range 
of wind directions and therefore reduce the integrated measuring time or 
decrease the random (category A) uncertainty. 

The benefits of an operational laser anemometer will be 

Reduction in direct costs of performance tests. 

V 

Reduction of the uncertainty of the individual power curve as well as of 
the entire wind farm. 

Lesser financial costs following the lesser uncertainty of the power curve 
(investment in wind power plants becomes less risky when output is better 
predicted; therefore a lower interest rate is asked by investors). 
With a laseT~anemometer as standard equipment on (large) wind turbines the 
logistics of performance (and other) tests will be minimised 
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MAIN GROUP 


COMPONENT 


UERTAINTY 


INSTRUMENTAL (1) 


Cup anemometer 






Calibration 






Mounting effects 


1-2% 




Operational effects 


0-2% 




Power sensor 


0.5-1.5% 




Air density 


0.5% 




Wind direction 


0.5% 




Precipitation 


- 




Data acquisition 


- 
- 


METHOD (2) 


Energy flux vs. hub height wind 


1-10% 




speed (incl. Turbulence) 




SITE CALIBRATION (3) 


Complexity terrain 


2-4% 




Correction size 






Regression RMS 






Wind direction 






variation 






Regression offset 






Distance 




OBSTACLES (4) 


Upstream wind turbines 






Upstream trees, houses 


0-3% 


EXTRAPOLATION (5) 


Sample uncertainty 


Potentially 






large 


DATA AMOUNT (6) 


Variability, power 


0-1% 




Variability, site cal. 





(1) Uncertainties directly related to calibration of instruments and their 
use; Most measurements consist of several measurements in combination. 
5 (2) Refers to the applied method itself (hub height measurements of wind 
speed in front of rotor etc.); the potential pitfalls are to some extend 
dealt with in section 3.1. 

(3) The present method, IEC 1996, highly relies on methods to re-calculate 
the wind speed measured at met-tower to wind speed at the wind turbine; 

10 with the IEC method site calibration uncertainties are potentially severe 
though very difficult to quantify. 

(4) Analyses, Frandsen 1994, indicate that obstacles (of sizes of the wind 
turbine) as far away as 1 km may influence the accuracy of the performance 
test. 
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(5) When evaluating plants of a larger number of machines, present practice 
is to test one or two turbines, which leaves sample uncertainty. 

(6) For both power variability (power curve measurement) and variability 
when performing site calibration there are (measurable) uncertainties 

5 Table 3.1 Summary of uncertainties in power performance testing. Only 

uncertainties related to "amount of" data are quantifiable. The indicated 
uncertainty levels are indicative . 

The operation of a wind turbine can be divided in four modes depending on 
10 the wind speed: 

Below cut-in 
cut-in to rated 
rated to cut-out 
15 above cut-out. 

The behaviour of the wind turbine in the two power production modes depends 
on its type, whether it is stall-, pitch or pitch/variable speed 
controlled. Considered here are only pitch and combined pitch and variable 
20 speed controlled machines. Both these two kinds of wind turbines apply 
active control mechanisms for limiting the output power. The same 
mechanisms can be utilised for production optimisation i.e. operating the 
wind turbine at maximum Cp. 

25 The basic objective of a control system is to reach a given value of an 
output signal through manipulation of an input signal. There are two 
separate ways of controlling a given system: open and closed loop. In the 
open loop approach the input variable is assigned a value based on a priori 
knowledge of static characteristics of the controlled system e.g. to have 

30 500 kW from a given wind turbine the wind speed is measured to 15 m/s and 
then the pitch angle should be 5 deg. and the rotational speed should be 30 
rpm based on pre-calculated Cp curves. The actual output is not measured. 
In the closed loop approach the pitch angle and rotor speed are determined 
on the basis of measurement of the output power. This ensures that 

35 deviations from the assumed characteristics and disturbances are 
suppressed. 

The closed loop performance is limited by several factors. It is seen, in 
figure 4, that if the phase shift through the loop is 180° and the loop 
4 0 gain is above 1 the fed back signal is actually added to the reference 

signal and the resulting system is therefore unstable. Also in the figure 
it is shown how a disturbance, d, can be fed forward, through H 2 , in order 
compensate for the disturbance before it lias changed the output through G 2 
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and Gi. It is of course necessary to be able to measure the disturbance if 
it is to be fed forward. 

The output power from a wind turbine has to be limited in order to limit 
the loads on the various components. The output can be limited e.g. by 
applying variable pitch and/or variable speed. 

For pitch control the basic control scheme can illustrated as in figure 5. 
The objective of the control system is to keep the torque constant (it is a 
constant speed machine), in the current case to keep the torque at 110 kNm 
for all wind speeds where it is possible. 

The setting of the pitch angle is controlled by a controller, typically a 
PI controller with electrical power as input. The electrical power is used 
for several reasons. It is among other things the power that has to be 
limited and it is simple, robust and cheap to measure. From a control point 
of view the electrical power is a more problematic signal to use. The basic 
problem is that instant torque on the rotor is determined by the instant 
pitch angle and the instant wind speed as experienced by the rotor and the 
power the is used for the control loop is the electrical power that appears 
after the torque has been converted by the gearbox and the generator. The 
controller then has to calculate the pitch angle that compensates for the 
deviation in power production and the pitch system has to move to the 
desired pitch angle in order to limit the power. Unfortunately the very 
high sensitivity to wind speed changes and the phase lag resulting from the 
conversion of wind power electrical power makes it very difficult to obtain 
a power limiting performance that is adequate. 

One of the classical and best ways of improving a control system is to feed 
forward the signal that causes the disturbance in order to compensate for 
the effect it will cause before it actually causes it. For the pitch 
controlled wind turbine looked at here The desired signal to feed forward 
is the instant wind speed experienced by the rotor. This signal is very 
difficult to obtain on a standard wind turbine because the anemometer on 
top of the nacelle is not accurate enough due to the disturbance from the 
rotor itself and it is not practical to have an array of anemometers in all 
directions to get an upstream measurement of the wind speed. 

The active control scheme can also be utilised for energy production 
optimisation. The pitch angle or rotor speed is then controlled in order to 
operate at maximum Cp. Various ways of implementing the optimisation are 
possible but the simplest is to use the wind speed and then, based on 
steady state characteristics, determine the optimal pitch angle or 
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rotational speed. The principle for a variable speed machine can be 
illustrated as in figure 7. It is desired to operate the rotor at a 
constant tip speed ratio in order to have maximum efficiency. The tip speed 
ratio is the ratio between the wind speed experienced by the whole rotor 
5 and the rotational speed of the rotor. It is therefore very desirable to 
know the wind speed. 

If the wind speed is to be used for control purposes there are several 
considerations that has to be taken into account, it is clear that the 
10 better the wind is known the safer it can be fed forward. 

The wind speed that is needed for power limit control is the wind speed 
seen {or to be seen) by the whole rotor. This wind speed is very difficult 
to obtain, because the wind turbine itself disturbs the flow and the wind 
15 speed is measured at one point. 



On the current wind turbines the cup anemometer is mounted on the nacelle. 
The wind speed measured there can not be used as wind speed input for power 
limit control because the wind speed at the anemometer is disturbed by the 
20 rotor itself. Instead of using an anemometer on top of the nacelle a laser 
anemometer can be used to measure the wind speed at a point (or several 
points) up stream. 

The measured wind speed upstream can then be used to predict the wind speed 
25 that the rotor will see at some point later in time. The wind speed 

prediction will have to take the delay, the spatial filtering of the rotor 
and correlation between the wind speed measured at one point and the wind 
turbine's position. The estimate of the wind speed can be derived in two 
ways. It can be done in open loop as in figure 8 or it can be done in 
30 closed loop as in figure 9. 

The estimated wind speed can then be used in the control algorithm. All the 
models involved are dynamic. The wind turbine model should be as simple as 
possible only taking low frequency (~1 Hz) dynamics into account. 
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Since the final goal of this investigation is to assess, with blade pitch control in 
mind, the possibility of using the upwind air velocity, measured from the wind 
turbine, a number of micrometeorological questions have to be answered. If we need 
five second of time, say, to carry out a change in the pitch, then we must know how 
well the measured wind speed, perhaps ten rotor diameters upstream, represents the 
actual wind speed at the rotor five seconds later. 

In the following we will address the questions: 

1. Eddies of smaller extension than the linear dimensions of the rotor are 
uncorrelated over these dimensions and of no importance to the overall forcing. 
Consequently, we must emulate the area averaging of the wind field by the rotor 
as well as possible in the upwind detection. We assume that this is carried out 
by a line averaging in the wind direction. What should the length of this line 
then be? 

2. What is the correlation coefficient between the upstream wind speed and the wind 
speed at the rotor? 

3. How do we interpret the correlation coefficient to tell us how well the measured 
upstream wind speed matches the wind speed at the rotor. 

Let the spectral tensor component for the streamwise, unfiltered velocity component 
u°(x) be $J tt (fc). Then spectral the tensor component of the spatially filtered velocity 
component u(x) will be given by 



In general this relation is given by the integral over the entire physical space 



Utl 



(*). 



(4.1) 



where the actual form of 5(fc) depends on the relation between u(x) and u°(x). 




(4.2) 



V 



The volume transfer function E(k) is then given by 



E(k) = \h(k) 



2 



(4.3) 



where 




(4.4) 



v 
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To obtain the usual one-dimensional spectrum F u (k\) of u as a function of the 
streamwise component of the wave vector fc, we carry out the double integration over 
the other two components of k, viz. 

= r dk 2 r dk z ^ uu (fc) = r r ^sc*)*^*). ^ 

J—oo J —oo J -oo J -oo 

The corresponding unfiltered one-dimensional spectrum F^(k) is of course given by 

= l°° dk 2 dk 3 *° uu (k). (4.6) 

J— oo J — oo 

The one- dimensional transfer function is defined as 

"<*>=it M7) 

In the special case of isotropy we have 

*2„(*)=^^. (4-8) 

where E(k) is the energy spectrum and in this case there is a simple relation between 
F%(k) and E(k)> namely 

(49) 

In this case ue assume that /i(x) = /i*,(x) is given by 

hL^) = \ T ^-2&{x 2 )8{x 3 ), (4.10) 

where L, as indicated in Fig. 19, is the characteristic length of the one-dimensional 
spatial weighting. Substituting in (4.4), we get 

h L (k) = e-^ L . (4-11) 

The volume transfer function thus becomes 

E*(*) = e- 2 ^. (4.12) 

Since E L (k) is independent of k 2 and A*, (4.5), (4.6) and (4.7) imply that the 
one-dimensional transfer function is equal to 7it(k) , i.e. 

H L {k i )=E L (k) = e^ L . (4-13) 
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We now assume that the mean wind and, consequently the Xi -direction is perpendicular 
to the rotor of the wind turbine. The spatial weighting function h(x) = h R (x) is 
therefore 



h R (x) = <S(xi) 



.1 



0, s/xl+xi >R 



(4.14) 



where R is the radius of the rotor (see Fig. 19). 

We substitute (4.14) in (4.4) which because of the delta function reduces to the 
double integral 



(4.15) 



Rotor 



This integral can easily be carried out by using plane polar integration variables: 



M*) = / rdr I exp(ir{* 2 cos<p + k 3 simp}) dip 

ntt Jo Jo 



= 2 v . ,. y , (416) 

where J„(s) is the Bessel Function of the first kind and of order v. 
The volume transfer function Er(Ic) is according to (4.3) 



Substituting this into (4.5), we obtain the filtered one-dimensional spectrum F*(fci). 
In order to be able to carry out the two integrations we must assume a particular form 
of the spectral tensor component $° u (fc) and here we assume that the turbulent velocity 
field is isotropic so that we can use (4.8). The integration could be carried out over 
k 2 and * 3 , but because of the form of 7i R {k) in (4.17) and $l u (k) in (4.8) it is 
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advantageous to change to plane polar integration variables 













cos# 




> =K x < 


> 






sin$ 


k 4 







(4.18) 



Fu^i) = P dk 3 P dk 3 S R (k)*° uu (k) 

V— oo J—oo 



{.kl+kl)B? 



roo roo 

= / dk 2 I o 

J — oo J -oo 

9 f oo E(jh 2 + K 2 ) 



47T 



(*? + *2 + * 2 ) 2 



(4.19) 



We cannot carry out the final integration over K without making a specific assumption 
s about the energy spectrum E(k) and here we chose the simplest possible, namely that 



E(k) = ae 2 ' 3 ^' 3 , 0 < k < oo, 



(4.20) 



where a« 1.7 is the dimensionless Kolmogorov constant and e the rate of dissipation of 
specific kinetic energy. 

Equation (4.9) implies that 



i?(*) = ^(i*i) = ^«* 2/3 i*r 8/3 - 

Substituting (4.20) in (4.19), we get 



(4.21) 



P 

Jo (* 2 * 2 +s 2 ) 17/6 



ds 



= ^a£ 2 / 3 |*|- 5 / 3 x #,,(*), 



(4.22) 



where 
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where \F2(a\b,c\x) is the Generalized Hypergeometric Function. 

In the derivation we have assumed uniform weighting over the entire rotor plane. 
However, since the largest part of the external torque from the wind is located far 
from the axis, we have also calculated Hr(Ic) in the case when the spatial weighting 
over the rotor plane increases with the square of the distance from the rotor axis. 
Again, it was possible to derive an explicit expression like (4.23), albeit a lot more 
complicated. 

Figure 20 shows Hji{k) t given by (4.23) together with a simple exponential expression 
exp(-0.75A;jR). Also shown is H R (k) in the case when the uniform weighting over the 
rotor is replaced with a weighting which increased with the square of the distance 
from the rotor axis. 

Assuming that the exponential expression is accurate enough to characterize the area 
-filtering by the rotor, we conclude that the lidar line "averaging should have a 
characteristic length L of about 

L « 0.375#, (4.24) 



irrespective of the weighting over the rotor plane. 

The question about how well correlated the upstream longitudinal velocity component is 
to the same velocity component just in front of the rotor of the wind turbine will be 
discussed in terms of the so-called longitudinal velocity coherence. If Taylor's 
hypothesis were strictly valid this coherence would be exactly one and we would expect 
perfect correlation. However, Taylor's hypothesis is of course only approximately true 
and we must consequently move a little outside the realm of standard assumptions when 
we deal with the longitudinal coherence. 

From an experimental point of view, deviations from unity longitudinal coherence are 
difficult to determine. The reason is that, even when making a strong effort, the line 
between the two observations points will in general only by coincidence be 
sufficiently close to mean-wind direction during the observation period. This 
direction is known only after the measurements have been carried out and if there is 
just a small angle between the mean-wind vector and the separation line between the 
anemometers the measured coherence will be dominated by lateral coherence which falls 
off with frequency much faster than the longitudinal coherence. We will therefore 
first make an attempt to understand the longitudinal coherence and the corresponding 
correlation on basis of a simple, serai-quantitative model and then estimate the 
contamination from lateral coherence due to misalignment. 
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Consider the covariance 

Ru(X>t) = ((u(x,0 - U)u(x + X y t + r) - I/)), (4-25) 

where u(x,f) is the instantaneous wind-velocity component along the mean-wind 
direction at point x and time t, [/ = <u(z,f)> the mean-wind speed, X the separation 
vector and r the time lag. Then the cross -spectrum is 



X*(X,») = ±J™ ga Ru(X,T)e- i '"dT 



(4.26) 



and the coherence defined as 



coh u (X,u) = 



Xu(0,u>) 



(4.27) 



This definition is general and can be used irrespective of the orientation of the 
separation vector X. 

In the case when X is along the mean-wind direction, i.e. X = Xi x , Taylor's 
hypothesis implies 

R»{Xi 1 ,T) = R»(o,T-^) (4-28) 
and, when substituting this in (4.26), we obtain 

x »(Xiu») = ± V ) e " iUTrfT = (429) 

We see that there is only a phase difference between x u (A'ii,u>) and Xu(0,w). 
Consequently the coherence, defined by (4.27), must be one, irrespective of the 
frequency, in cases where we assume Taylor's hypothesis. 

Another way of stating this consequence of Taylor's hypothesis is to interpret (4.28) 
and note that u((x + X)ti , t) is perfectly correlated with u{xiut-X/U)= the upstream, 
streamwise velocity component at a time so much earlier as it takes a disturbance to 
travel the distance X with the speed U . In other words, 

({i,(«<„«-^)-t/}{«(l» + xji..«)-v}) = n( x *4) 

= R v [0,0) = o 2 u , ( 43 °) 

where at is the variance of u. 
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We will modify (4.29) to account for the decorrelation due to the turbulence by the 
generalization (Kristensen, 1979) 



(4.31) 



where the exponential 'decay factor* is a decreasing function of A'. The 'eddy 
mean- lifetime* T(u>) is assumed to be a function of the eddy size A ~ 2itU/u)* . 



Mann (1994) found that 



TM-^(^n-(i.?4-[^n} 



-1/2 



(4.32) 



where T«3 is a dimensionless constant, k « 0.4 the von Karman constant, u, the 
friction velocity and z the height, brings his spectral model in good agreement with 
experimental data over land_and_ov.er water. 

for very numerically very small and very large values of w/U we have 



3B 



1.2 



(*)". 



M -» 00 



(4.33) 



-2/3 



where 



B(a,y) = 



r(x) xTjy) 



T(x + y) 

is the (Complete) Beta Function. 

As we are concerned with rather large eddies, we shall use the asymptote 

. -l 



(4-34) 



7TM«L5i(g.)- 



(4.35) 



All we now have to do to determine the correlation coefficient with longitudinal 
separation 

piXti) - — a — 



(4.36) 



is to apply the transform 

Ru(Xi u r)= f°° Xu(Xi x% u)e>" T <b> y 

J —CO 



(4.37) 



♦ Actually we cheat a little here and in the following and use Taylor's hypothesis to translate eddy size or reciprocal 
wave number to frequency u/. 
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with r = X/U 9 to Xu(Xi u tj) given by (4.31): 



«■(**•*) - £*(-!t7{W)«*->*- 

where F U (A:) is defined by (4.5). 

We assume, according (see Fig. 20) that 

F u (k) = e-°' 76kR F^(k) (4.39) 

and we will use the von Karman spectrum with a length scale equal to z, consistent 
with Mann (1994) , 77e~ 

a 2 z 

^ ik) = 5(1 A 1/3) (1 + (fo)2)»/6 ' (4 - 40 > 
Further, we will assume that the rotor radius R is half the hub height z. 
Substituting in (4.38), we get 

where 



- 3^F 2 (l;^ ; -{|} 2 )-^(|) 1/3 Y. 1/3 (a (4.42) 

where Y„(s) is the Bessel Function of the second kind and of order v. 
We can now determine the correlation coefficient p(Xii): 

The correlation function is shown in Fig. 21. 

If the normal to the rotor plane, the direction of the yaw, is not parallel to the 
mean wind direction, a gust detected at the distance X will not travel along the yaw 
direction and will therefore not hit the center of the rotor. The correlation will 
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consequently not be perfect between the upstream wind velocity and the velocity at the 
rotor at the time X/U later. Let the angle between the yaw direction and the 
mean-wind direction be a. Then, assuming Taylor's hypothesis, the correlation in 
question is that between velocity components laterally separated the distance 
Xtan(a)«Xa. By simple modelling we can estimate this correlation and we will 
consider both point velocities and velocities which are averaged over the rotor plane. 
We start with point velocities. 

The spectral tensor component for the longitudinal velocity is given by (4.8) and here 
we use the energy spectrum which can be derived from the von Karman spectrum and the 
relation (4.9) : 

Elk) - I °* z {k2) * , . (4.44) 

The spectral -tensor component can be written 

.« fn 1 °V (445) 
*" W " 6ir 5(3/2,4/3) {1+ (*? + Jbf + &f)z 2 }"/« ' 

The covariance of the two unfiltered velocity components with the lateral separation 
Xa becomes 



6tt 5(3/2,4/3) /_„ J 0 Jo {1 + (*? + K 2 U 2 }"' 6 

~ 3 5(3/2,4/3) /„ (1 + .*+«»)»/• /o J °U J 

- ^r^r)o^ rf5= ^° M ' (446) 

where p°(Xai 2 ) is the the correlation function, given by 



where K u {s) is the a Modified Bessel Function (convergent for s oo) of order i/. 

Here and in the following we consider the typical situation where the hub height is 
twice the radius of the rotor, i.e. z = 271. The last result can therefore be written 
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\X«i 2 ) = |^{f(473j(S) Kl/3 (S) 

When the velocity signal is averaged over the rotor we must replace *° uu (k) by 

, M o5 2 "2* g(vg^) (4.49) 

*«.(*) = 2 3^5(3/2,4/3) {1 + 4 (fc? + fc| + fc 3 2 )fl 2 } ,7/6 - 

The covariance becomes 



o ^2p3 /-oo /-oo z-2* j2(/CA'a)exp(i/'!:A'acos^) 

^hmbmL^L KdK L » (i + 4(*?+«')*>}"/« 



J?(5)J 0 ([Xa/fi)3) ^„ (4.50) 
{l + 4s 2 } 7/3 



The correlation coefficient 

/■« J?(s)J 0 ([Aa/.R]s) / r j|(£) (4.51) 
p(Xat 2 ) = J o {1 + 4s2} 7/3 sds / y 0 {1+45^ 

can be evaluated by numerical integration as a function of Xa/R. Figure 22 shows the 
correlation when the velocity is unfiltered and when it is area filtered. 
Until now we have only discussed the correlation coefficient po between the upstream 
velocity and the velocity at the wind turbine. It seems tempting to assume that if Po 
is close to one then a large turbulent excursion of the upstream velocity will give 
rise to an excursion which is similar in shape and magnitude at the wind turbine at a 
time X/U later. How can this assumption be quantified? 

Let us call the upstream velocity u, and the velocity at the turbine u. Let the 
variance *l and the correlation coefficient p 0 be given, and let us assume that u, and 
u have a joint Gaussian probability density 

^(-27T^ {M '- 2po " ltt + M2} ) 



(4.52) 
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We may recast this expression as a product of the probability density function Po(«i) 
for ui and the conditional probability density Pi(ui;u) for u. We are particular 
interested in Pi which becomes 

We see from this expression that the conditional mean value of u is PqU\ and not u x . 
The implication is that if p 0 < 1, the mean (and the most probable) value of the 
velocity u at the wind turbine is less than ui if tii is positive. A given upstream 
excursion will be suppressed on its way to the wind turbine. 

Since u is biased in such a way that its most probable value is uq = p 0 u\ , it seems 
reasonable to use this value for controlling the rotor. Consequently, we may ask the 
question: what is the probability rpi that |u - uo|/|ito| is less than a small positive 
- number - e? 

This probability will depend on p 0 > £ aui d t ^e excursion in units of the standard 
deviation uo/a u . 

uo+<|u 0 | , 1 l\ 

uo-<|uo| 

where 

erf(x) = f'e'^ds (4.55) 
V* Jo 

is the Error Function, shown in Fig. 23. 

When x > 2, the Error Function can conveniently be approximated by 

erf(x) « 1- ^— . (4.56) 

We might get a better estimate of the excursion u at the wind turbine if we measured 
the wind speed at several places upstream (a range-gating system) . Let us assume that 
we measure the wind speeds at two upstream positions and detect the excursions u\ and 
u 2 . Let the correlation coefficients between u and u x and tt 2 be pi and p2 • 
respectively, and let the correlation coefficient between u\ and u 2 be p. Further, let 
us assume that u, u\ and u 2 have a joint Gaussian probability density. Then the 
conditional probability for u, when u\ and u 2 are given, becomes 



1 / 1 - P 2 1 
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We see that in this case the most probable value of u becomes 

Pi ~ PP2 , P2 -pPi .„ (4.58) 

The probability V>2 that |u - uol is less than c|uo| is in this case given by 

J) = erf (\/ 2(1 + 2ppip2 -ti-pl)^) " (4 * 59) 

In the limit when p = 1 and p, = p 2 2 po we obtain the same result as in the case when 
there is only one upstream measurement with t*o = Pofai +«a)/2- 
_.It might be interesting to see if we get a better„estimate with two upstream 
measurements if pi = p 2 = Po and p < 1 . In this case (4.59) becomes 

=erf (/^^^ € ^) ' (460) 

Since 

1+P > 1 (4.61) 
2(l + p-2pg) " 2(1 -pg) 

we find 02(po, Po»P> £*u 0 /a u ) > tMpoiC u 0 MJ • 

If we also let p be very close to unity it might be useful to evaluate the difference 
^ = ^(po,P,e,W^)=^(Po.Po,P^,WO-^i(po^,«oK) to the first order in 1 - p. 
The definition (4.55) of the Error Function implies 



Sxp — erf 



= erf 



erf 



( I 1+P C M \ *Tf( . 6 „ K|N j 

^/2(l + p-2pS) <7,J V^l-pg) ». / 

f (i-P)pg i £ N) 

\4(1-pS)v/2(1-pS) <W 



1 (l-p)pg Kj (4.62) 
" 2v/2lF (1 - PS) 3/2 ' 

As an example, let p„ = 0.9, < = 0.2 and |uo|K = 2. Then (4.54) and Fig. 23 shows that 
0,(po.e,uo/*»)«O.64. From (4.62) we get the increase t* « 0.78(1 - p) when two upstream 
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velocities are used. We expect that the spatial configuration implies that p<p 0 . 
This means that the probability increase is less than about 0.08, corresponding to a 
12*/, 1 ' improvement 9 ' . 

Several different systems for localised measurement of wind velocity have 
been devised. Most practical applications are based on the Doppler effect 
with optical hetero- or homodyning. 

A laser is not absolutely mandatory, but it is very difficult to get good 
measurements with an incoherent light source. Thus only lasers are 
considered for the light (radiation) source. The lasers considered for the 
present application are listed in Table 1. 



Table 1. Potential and actual lasers for long-range laser anemometry. It is 
noted that the holmium-l-aser is attractive for wind measurements; however; — 
the price and complexity of current lasers prohibit the application in 
low-cost systems. Future developments may change this. For the reference 
beam configuration it is mandatory that the intrinsic laser noise (not the 
photon noise) is extremely low. 



Laser 



Wave- Power Coherence Efficien- Compact Eye- safe Price 

length output proper- cy and 
urn ties and robust 
noise 



C02 



Tin, Ho : YAG 



Nd: YAG 



GaAs 



Argon 



10.6 



2.1 



High 



High 



1.06 Medium 

(0.53) high 

0.85 low 

0.488/0.5 medium 
14 



Good > 10 % yes yes 

(waveguid 
e) 

medium - 1 % {not yes 

currently 
) 

medium - 1 % partly no 



medium - > 10 % yes 
low (very) 



good 



no 



0.01 % no 



medium - 
low 

very high 

medium 

low 

high 



At the current state of technology only the C02 laser is a viable option 
for an efficient long-range system with a potential for being manufactured 
at a reasonable cost. 



Solid-state lasers may become attractive in the future. An applicable 
solid-state laser should have a wavelength in a transmission window of the 
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atmosphere and should preferably be eye-safe. Solid-state lasers that 
fulfil these requirements are based on YAG (yttrium-aluminum-garnet) as 
host material and are doped with a material defining the optical 
wavelength. This may be rare earth materials. The price of such lasers 
5 (Tm, Ho r YAG and equivalent) is currently very high (~ USD 10,000 and up) . 
NdtYAG lasers pumped by semiconductor lasers could be a viable alternative 
because they are more likely to be available at a reasonable price. 
However, they are not eye-safe. 

10 GaAs semiconductor lasers can be produced at very low prices. Advances 

within arrays of phase locked lasers may even make it possible to produce 
lasers with adequate output power generated with a very high wall-plug 
efficiency (considerably larger than 10 %) . It is, however, not likely that 
low-cost eye-safe semiconductor lasers will be available in the near 

15 future. 



The argon laser has been the "workhorse" of short-range laser doppler 
anemometers. However, it is not eye-safe and the wall plug efficiency is 
very low. The wavelengths that can be obtained are generally too short for 
20 the present application. 

The optics serves to reshape the beam from the laser, to transmit it to the 
measuring region, to collect the scattered light and to direct it to the 
detector, possibly with a properly mode-matched reference beam. The optical 
25 components can be classified as either refractive/reflective or 
dif f ractive . 

Refractive (and reflective) : change of direction of propagation upon 
passage from one medium to another. Only simple optical operations are 
30 possible in each element; the mechanical construction supporting the optics 
is in general expensive; the volumes of the elements are as a rule large. 

Diffractive (incl. holographic): surface or volume structure causes 
interference in the form of diffraction that may mimic conventional optical 
35 components - but it is also possible to implement components with no 
refractive equivalent; multiplexing is possible. These elements are 
expensive as single elements, but (potentially) of low cost in mass 
production. 

40 The type of optics very much depends on the optical wavelength. Glasses 
that are transparent in the visible and near infrared cannot be used at 
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longer wavelengths. ZnSe is now the preferred material for refractive 
optics at the long infrared wavelengths. 

Optical components must usually be fabricated to tolerances much smaller 
5 than the optical wavelength. Intrinsically this makes it cheaper to produce 
infrared optics (often counterbalanced by the more expensive materials that 
are also difficult to process) . Reflective metal optics is attractive for 
especially C0 2 laser applications. The reflective optics can be combined 
with diffractive structures made as relief patterns on the metal surfaces. 
10 Beamsplitters and isolators are delicate and expensive for the wavelength 
of the C0 2 laser. 

Scanning a plane (or strictly speaking: part of a sphere) may be performed 
with sets of mirrors. This is sometimes done to get a spatial average or to 

15 obtain different velocity components of the mean wind. The principles are 
simple but the mechanical layout of existing implementations is bulky. 
Different ranges can in principle be covered by incorporating a zoom 
facility in the optics. If scanning is to be incorporated, the beam 
expanding optics may be designed to incorporate this facility. This would 

20 provide a much more compact and robust implementation, but could not be 
implemented in a compact form with off-the-shelf components. Ray-tracing 
calculations indicates that with diffraction limited performance and only 
two optical elements it may be very difficult to obtain more a scan larger 
than 10 % of the mean range or angular position. 

25 

In order to obtain information about the movement of a given volume element 
in a fluid, it is necessary that the propagation of light is perturbed in 
the volume. Ideally the light should propagate unperturbed to the measuring 
volume. In the measuring volume light should be redirected to the receiver. 
30 The situation is rarely that ideal: the propagation of light is affected by 
refractive index variations. These may be generated by turbulence causing 
temperature fluctuations that again give rise to small fluctuations of the 
refractive index of the atmosphere. For shorter ranges (< 100 m) in the 
open air this effect is normally negligible. 

35 

Scattering by turbulence induced refractive index fluctuations cannot 
generate the desired return signal. The spatial distribution of the 
refractive index fluctuations is identical to the distribution of the 
turbulence induced temperature fluctuations. The spatial scales are larger 
40 than the Kolmogorov or inner turbulence scales, which in the atmosphere are 
larger than a millimeter. The scattered light is essentially confined to an 
angular region given by the ratio of the optical wavelength to the size of 
the scattering structure (weak interaction assumed), i.e. 
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a m m^lK mt (17) 

which implies that light is only scattered within a narrow cone in the 
5 forward direction. 

However, there are in general many small particles in the atmosphere. They 
occupy a large size range from below 10-2 urn to above 10 urn. The size 
distribution is often described by the so-called Junge particle size 

10 distribution. In a standard atmosphere the background aerosol concentration 
is about 400 particles/cm 3 . The concentration close to the ground is 
usually much higher especially the concentration of large particles. The 
concentration and distribution to a large extent depend on the weather 
conditions and emissions from the ground. 

_15 - - 

Light scattered by small particles has been the subject of research for 
more than a hundred years. A well established theory exists for spherical 
dielectric particles. It is essentially so that for particles much smaller 
than the optical wavelength, backscattering and f orwardscattering are about 

20 equal. As the particle size increases, forwardscattering increases much 
more than backscattering. (see, e.g., Bohren and Huffman.) 

In the reference beam backscattering system the transmitted and the 
received light appears in the same spatial mode (except for the direction 
2 5 of propagation) . This fact has the consequence that for a focused 

continuous wave system the collected light power is independent of range 
and transmitter/receiver aperture as can be seen below. 

The collected light power is given by the following simple expression: 

30 P r ={a)cV%-^r, < 18 > 

A 0 z 

where (a) is the mean scattering cross-section, c the particle 
concentration, V the volume of the measuring volume, P t the transmitted 
laser power, Ao the cross-section of the measuring volume, A c the 
receiver/transmitter area and z the distance to the centre of the measuring 
35 volume. The parameters for a Gaussian beam are as follows: 

(19) 

n ( 

Az = r Q 2 k - —X 



4 



z 
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from which Ao, A c , and V are evaluated and substituted into eq. 18. This 
gives that the received power is given by 

P r =^- y k = 2nlX, (20) 
k 

i.e. proportional to the following product 

5 

mean scattering cross-section x particle concentration xwavelength. 

The backscattering of light is calculated in a single spatial mode 
expressed in terms of photons/sec assuming the scattering volume is 

10 illuminated with a beam that propagates in the same spatial mode - except 
for the direction. This appears to be very relevant when the power 
requirements for a reference beam configuration is investigated. The 
results as a function of the optical wavelength are shown in figure 11. It 
is noted that the curve peaks at around 10 urn almost irrespective of the 

15 specific particle size distribution. The simple C02 system has its 

wavelength at 10.6 \m\ A constant refractive index is assumed. This is 
probably a good approximation for a background aerosol that mainly consists 
of various metal oxides. For water particles the refractive index increases 
drastically with increasing wavelength, which enhances the backscattering 

20 cross-section relative to the values obtained for the visible part of the 
spectrum. The backscattering may be 2 - 3 times higher at a humidity close 
to 100 % relative to 0 % humidity. 

Most particles in the atmosphere are not spherical. The exact shape has a 
25 small impact on the light scattering characteristics if the particles are 
smaller than the optical wavelength. The analysis of light scattering from 
large particles of arbitrary shapes is in general difficult - or impossible 
- to calculate. However, there is good experimental evidence that average 
values for light scattering can usually be estimated quite well assuming 
30 spherical particles. 

Light must propagate through the atmosphere from the transmitter to the 
measuring volume; and the scattered light must propagate back to the 
receiver. Ideally light would propagate completly unaffected by the 
35 atmosphere - except for the light scattering in the measuring volume. The 
conditions are actually generally quite close to this situation (this is 
essential for the success of laser anemometry in general) . However, at 
longer ranges and in adverse weather conditions three phenomena that may 
have an impact on the propagation by: 

40 
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Extinction caused by absorption and especially by scattering. In clear 
weather the effect of this is marginal, but may be substantial in heavy fog 
or rain. There is a small absorption of light by C0 2 in the atmosphere at 
the wavelength of the C0 2 laser. For very long-range applications 
5 (kilometres) an isotope of C0 2 has been used for the laser to obtain a 
slightly de-tuned wavelength. For shorter ranges (such as in the present 
case) the absorption is negligible. The extinction caused by heavy fog may 
from 15 - 60 dB/km at 10.6 urn. At the ranges of interest in the present 
context this would only correspond to a variation in extinction by a factor 
10 of 2 - 3. This would partly be compensated by increased scattering under 
such weather conditions . The effects are considerably larger at shorter 
wavelengths . 

Perturbation of the wavefronts by turbulence induced refractive index 
15 variations in the atmosphere will be negligible at a wavelength of 10.6 urn 
and at the ranges considered here, but is considerably larger at shorter 
wavelength. Data on the wavefront perturbations caused by strong rain and 
fog are not available but based on heuristic arguments it is anticipated 
that the effect is small (i.e. does not prevent useful measurements) . 

20 

The photodetector serves to convert the optical signal to an electronic 
signal. Phototdetectors relevant to laser anemometry are all of a type 
where photons are absorbed whereby free charge, carriers are produced. The 
carriers (electrons) may be released from the material in which the photons 
25 are absorbed generating a current in an external circuit, or they may 

generate electron-hole pairs in an pn-junction also causing a photocurrent 
in an external circuit (photovoltaic effect) . The free carriers may change 
the conductivity of the responsive element, which is exploited in 
photoconductive detectors. 

30 

Before the characteristics of some detectors is considered, it is 
appropriate to consider the noise, i.e. random fluctuations of the signal 
that limit the performance of the system. The different types of noise that 
may play a role in the present context can be classified as follows: 

35 

• Light is detected in quanta with a random distribution (intrinsically a 
Poisson distribution) . 



• The detection process itself is a random process: not all photons 
40 impinging on the detector will generate free carriers and it is random 

which photons actually do generate charge carriers. 
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• Even without illumination free carriers will often appear and they 
generate noise. 

• The random generation and recombination is a random process. 

5 

• Thermal noise caused by the random motion of free charge carriers in the 
responsive element and/or in the subsequent electronics . 

• Lattice G-R noise; the random generation and recombination of free 

10 charge carriers may cause noise especially in photoconductors used at 

infrared wavelengths. 



15 



• 1/f-noise. The physical mechanism of this type of noise notably found in 
solid-state detectors is not well understood. 

Photodetectors for the detection of modulated light can be classified into 
the following categories: 



20 



25 



• photomultiplier tubes; 

• p-n photodiodes; 

• photoresistors ; 

• % drag' detectors; 

• thermal detectors; 

• phyroelectric detectors. 

Various types of preferred detectors are compared in Table 2. 
Table 2. Comparing different types of photodetectors. 





p-n diode 


avalanche 
diode 


photo res is tor 


PMT 


Quantum eff . 


High 


high - medium 


high 


low - medium 


Dark current 


Medium 


medium 


medium - high 


low 


Anpl ifi ca ti on 


1 


medium 


medium 


very high 


Speed 


Fast 


fast 


fast 


fast 


Spectral 


Very broad 


Broad 


very broad 


narrow 


region 


(different 




(different 


(visible + 




materials) 




materials ) 


adjacent) 



30 For the continuous signal the processor should not be gated but should be 
of a sliding window type (e.g. a frequency- or phase locked loop, or a 
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spectral analyser that operates in real time) . For the present application 
it is desirable that the processing incorporates a velocity estimation 
based on both the current measurement and the history the past evolution of 
the velocity (i.e. the likelihood of a large difference between two 
5 consecutive measurements is much smaller than the likelihood of a small 

difference - assuming the temporal resolution is smaller than the scale of 
the velocity fluctuations). A so-called Maximum A Posteriori (MAP) 
procedure is an integral part of tracking processors. A Kalman type of 
filtering of the temporal sequence obtained with a spectrum analyzer 
10 processor is also feasible. 

An outline of a laser anemometer for control of wind turbines is given in 

this section. 

On the basis of the discussion in the previous section it has been decided 
15 to consider a backscattering reference beam system based on a cw C02 laser. 
The system shown inTf Igure 12 . However, the price of a system" based on off- 
the-shelf components would be rather. Thus dedicated low- cost components 
would be needed for a system for actual wind turbine control. 

20 Considerable cost reductions are possible by introducing the following 
elements : 

• a custom designed waveguide laser; 

. combined reflective and diffractive metal optics, which could be 
manufactured on a precision turning machine. Polarisation elements, 
25 necessary for a reasonably efficient use of the laser power, could be 

incorporated by deep surface relief structures. 

A very simple system - which has received surprisingly little attention in 
the literature - could be based on a so-called autodyne configuration: the 

30 collected light power is fed directly into the laser itself and will 

therefore perturb the power output. The perturbation is monitored at, e.g., 
the rear end of the laser (figure 13). It is even conceivable to implement 
the laser as a tapered waveguide with an output window that serves both as 
laser mirror and focusing element. The tapering is done to reduce size and 

35 enhance laser gain. A high internal gain is important for the autodyne 
configuration seen in figure 14. 

The optics may be implemented in a form that allows for adjustment of the 
location of the measuring volume - possibly a conical scan or measurements 
40 at different ranges. Reflective optics would be preferable. The beam 
expander could be implemented as indicated in figure 15. Ray-tracing 
calculations indicate that diffraction limited operation with only two 
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elements may be very difficult if scanning capabilities must be 
incorporated. A combined reflective and diffractive system appears to have 
a better performance. 

A window will be necessary to protect the optics. This could be made of 
either germanium or zinc selenide. 

A likely signal processing scheme could be based on a frequency locked loop 
(described in section 5.1). A block diagram is illustrated in figure 16. 

The primary specifications are mean measuring distance and spatial 
averaging length. On the basis of the discussion in section 3 and assuming 
a rotor diameter of 60 m, the following values are selected: 



15 Mean measuring distance, L 150.0 m 

Spatial averaging, 2Az: 22.5 m 



This gives the following values: 
20 Focal diameter, dO: 17.4 mm 

Transmitter diameter, d t : 116 mm 

as illustrated in figure 17. 

25 The calibration of the system is given by (frequency/ velocity) : 

/ C rt , MHz 

^ = — = 0.189 — — . (21) 
v 2/r m/s 

A simulated detector signal is shown in figure 18. 

30 

There are a number of causes of the uncertainty of the measurement: 

• Intrinsic signal fluctuations associated with the random distribution of 
particles and/or turbulence. 

• Intrinsic signal fluctuations associated with turbulence on spatial 
35 scales smaller than the size of the measuring volume or on spatial 

scales smaller than the averaging time. 

• Photon noise 

• Thermal noise 

• Intrinsic fluctuations of the laser (laser noise) 



SUBSTITUTE SHEET (RULE 26) 



WO 98/42980 PCT/DK98/00125 

44 

• Perturbation of the beam propagation caused by random fluctuations of 
the refractive index of the atmosphere. 

The statistical uncertainty caused by the intrinsic signal fluctuations has 
5 been investigated by a number of researchers (originally by George and 
Lumley) . For the present application, the relative variance is given by 

var{v,} _ varfro,} a 1 

T~ ~ 2 ~ K A V ' 

v r0 a> d Ao)T <» 

where k is the relative spectral broadening. In the present case this 
broadening is given by the relative turbulence level on scales smaller than 
10 the length of the measuring volume. The minimum relevant averaging time is 
the transit time through the measuring volume, i.e. 

T = 2Az/v 0 . (23) 

Substituting eq. (23) into eq. (22) and utilising the fact that the 
relative spectral broadening is equal to the relative turbulence level 
15 yields 

var{v, } _ k 



(24) 

v., 



JO 



2 2AzC 



For the specific example, the relative uncertainty is 0.16 %. This must be 
considered to be the achievable minimum uncertainty. 

The condition for photon noise limited operation is established first. This 
implies that the photon noise dominates over the thermal noise. The thermal 
noise is given by 



« (v tn 2 ) = 4k B TRB f 



(25) 



where kB is Boltzmann' s constant, T the absolute temperature, R the 
resistance and B the bandwidth. The photon noise would generate a noise 
signal given by the following expression 

30 

(v pn 2 ) = £ 2 hvPB (26) 



where Q is the detector sensitivity (volts/watt), hv the photon energy, P 
the light power and B the bandwidth. (The expression is obtained by 
35 converting simple Poisson photon statistics to equivalent power 
fluctuations . ) Dominating photon noise implies that 
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C 2 hv 



For realistic values of resistance (« 100 ohm) and sensitivity ( > 10 
volts/watt) the inequality is fulfilled for power in the milliwatt range, 
5 which is no serious requirement to the laser and still low enough to stay 
within the maximum power level of newer non-cooled photoconductors . With 
the inequality fulfilled the signal-to-noise ratio is exclusively given by 
the collected light power and no detection scheme can provide a better 
signal-to-noise ratio (the only free parameter is the quantum efficiency, 
10 which for good detectors can only vary with a factor of 2). 

The following quantities are assumed to be given: 

• Five watt of laser power from a waveguide laser (a waveguide laser gives 
a-more compact and robust implementation) . 

15 • A negligible laser noise level. (In reference beam detection it is often 
the laser noise that limits the performance. However, it appears that 
single- frequency C0 2 lasers used in coherent lidar systems have such a 
low noise level in the spectral region of interest that the noise of the 
laser is negligible. This is consistent with experimental observations 

20 at Riso in connection with plasma diagnostics experiments using 

reference beam detection. 

• The collected light power is 10" 10 W. 

• An overall quantum efficiency of 40 %. 

• A detector with a detectivity D* of 1000 cmHz 1/2 W" 1 and an area A of 1 
25 mm*" . 

• A detection bandwidth B of 5 MHz. 



The signal-to-noise ratio is then found to be larger than 20 dB, which 
implies that the photon noise will have a negligible impact on the 
30 statistical uncertainty of the measurement. 

A system for on-line control could either be mounted on the top of the 
generator/gearbox casing or it could be mounted in the hub of the rotor. In 
the first case it is most likely necessary to incorporate a shutter 
35 facility that block the laser beam while a blade intersects the path of the 
beam. This facility serves to eliminate the very strong illumination of the 
detector by light reflected off the blades of the turbine. A simple 
rotating shutter synchronised with the main rotor could do this. 
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Mounting the laser anemometer in the hub would eliminate the blade 
intersection problem. However, power and signal connections would be more 
complicated. Some new large wind turbines do have a hollow main shaft in 
which electric wires can be placed or the laser beam could be transmitted 
5 through the hollow axle. A conical scan may be achievable by tilting the 
optical axis slightly with respect to the axis of the main shaft. 
Information about wind shear may be obtained in this way; thus making it 
possible to control the pitch angle in relation to the angular position of 
a given rotor blade. 

10 

The applicability of laser anemometers in conjunction with wind turbines is 
investigated. The main conclusions of the initial investigations are: 
Relevant measurements for both power calibration and control can be 
performed. 

15 

A reference beam C0 2 laser system is preferable to" the common LDA based on 
visible or near infrared light. 

A dedicated laser system is needed for a compact robust instrument. 
20 An autodyne system is very attractive in terms of robustness and 

(potential) low cost; however, the performance in relation to systems with 
external reference beam detection needs to be investigated. 

The effect of sever weather (rain, snow, strong optical turbulence etc.) 
25 need further investigation. 

In relation to the relevance of such systems it is noted that the industry 
estimates a 12 % annual growth rate for wind turbines that would benefit 
from a good feed forward control system and the number of wind turbines 
30 will be about 3000 units/yr. five years from now. 
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CLAIMS 

1. A wind turbine comprising a wind velocity measurement system for 
determination of air velocities in front of the wind turbine and having an 
output for provision of determined air velocity values. 

2. A wind turbine according to claim 1, wherein the wind velocity 
measurement system comprises a laser anemometry system having 

a laser for emission of a laser light beam for illumination of particles in 
the air, 

a detector for detection of light emitted from particles that have 
interacted with the illuminating laser light, and 

a signal processor connected - to the detector for generating a velocity 
signal corresponding to the velocity of the particles based on the detected 
light. 

3. A wind turbine according to claim 2, wherein the laser anemometry system 
is a Michelson interferometer type laser Doppler anemometer. 

4. A wind turbine according to claim 2, wherein the light emitted from 
particles in the air is transmitted to the laser for perturbation of the 
power output of the laser and the detector detects the perturbed power 
output. 

5. A wind turbine according to claim 4, wherein the laser comprises a 
tapered waveguide with an output window constituting a laser mirror and a 
focusing element. 

6. A wind turbine according to any of claims 2-5, wherein the signal 
processor is a frequency locked loop. 

7. A wind turbine according to any of the preceding claims, comprising 
adjustable pitch blades, 

a pitch controller for setting blade pitch and having an input, the blade 
pitch being set according to a signal provided at the input, and 

a control system for controlling the output power generated by the wind 
turbine by control of blade pitch and having 
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a control output that is connected to the input of the pitch 
controller/ 

5 a velocity input that is connected to the output of the velocity 

measurement system, and 

a pitch processor for calculation of the desired blade pitch based on 
determined air velocity values and generation of a corresponding pitch 
10 control signal that is provided at the control output, 

the determination of air velocities in front of the wind turbine 
compensating for delay in blade pitch setting so that the actual blade 
pitch setting corresponds to the actual wind velocity at the blades during 
15 operation of the wind turbine. 



8. A wind turbine according to claim 7, wherein the pitch processor is 
adapted to generate a pitch signal so that the tip speed ratio is 
substantially equal to a predetermined ratio. 

20 

9. A wind turbine according to claim 7 or 8, wherein the pitch processor is 
adapted to generate a pitch control signal so that the tip speed can not 
exceed a predetermined maximum velocity value. 

25 10. A wind turbine according to any of claims 7-9, wherein the pitch 
processor further comprises a data storage for storage of a first 
predetermined function of shaft torque as a function of blade pitch and 
wind speed for the wind turbine and wherein the pitch processor is adapted 
to calculate the desired pitch utilizing the first predetermined function. 

30 

11. A wind turbine according to any of claims 7-9, wherein the pitch 
processor further comprises a data storage for storage of a second 
predetermined function of normalized power output C p as a function of blade 
pitch and tip speed ratio for the wind turbine and wherein the pitch 

35 processor is adapted to calculate the desired pitch utilizing the second 
predetermined function. 

12. A wind turbine according to any of claims 11, wherein the pitch 
processor is further adapted to determine the pitch that for a given tip 

40 speed ratio substantially maximizes the normalized power output C p . 

13. A wind turbine according to any of claims 1-6, comprising 
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a control system for controlling the output power generated by the wind 
turbine by control of a load that is connected to the power output of the 
wind turbine and having 

5 a control output that is connected to the load, 

a velocity input that is connected to the output of the velocity 
measurement system, and 

10 a processor for calculation of the desired load value based on 

determined air velocity values and generation of a corresponding load 
control signal that is provided at the control output, 

the determination of air velocities in front of the wind turbine 
15 compensating for delay in load value setting so that the actual load value 
setting " corresponds to the actual wind velocity at~theTblades during 
operation of the wind turbine. 

14. A system for determination of the power curve of a wind turbine 

20 according to any of claims 1-6, comprising a power measurement system for 
determination of the output power generated by the wind turbine. 

15. A method of wind turbine power curve determination comprising the steps 
of 

25 

mounting on the wind turbine a wind velocity measurement system for 
determination of air velocities in front of the wind turbine, 

connecting a power measurement system for determination of the output power 
30 generated by the wind turbine between a power generator of the wind turbine 
and a load receiving the power generated by the wind turbine, and 

recording corresponding values of wind velocity from the wind velocity 
measurement system and output power from the power measurement system. 

35 

16. A method of controlling output power from a wind turbine, comprising 
the steps of mounting on the wind turbine a wind velocity measurement 
system for determination of air velocities in front of the wind turbine and 
providing determined velocity values to a power control system of the wind 

40 turbine. 
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